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Abstract: We present the calculation of the dominant next to leading order QCD correc- 
tions to Higgs boson production in association with three jets via vector boson fusion in 
the form of a NLO parton-level Monte Carlo program. QCD corrections to integrated cross 
sections are modest, while the shapes of some kinematical distributions change appreciably 
at NLO. Scale uncertainties are shown to be reduced at NLO for the total cross section 
and for distributions. We consider a central jet veto at the LHC and analyze the veto 
probability for typical vector boson fusion cuts. Scale uncertainties of the veto probability 
are sufficiently small at NLO for precise Higgs coupling measurements at the LHC. 
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1. Introduction 



One of the primary goals of the CERN Large Hadron CoUider (LHC) is the discovery 
of the Higgs boson and a thorough investigation of the mechanism of electroweak (EW) 
symmetry breaking |2[. In this context, vector-boson fusion (VBF) has emerged as a 
particularly interesting class of processes. Higgs boson production in VBF, i.e. the EW 
reaction qq QqH, where the Higgs decay products are detected in association with 
two tagging jets, offers a promising discovery channel @ and, once its existence has been 
verified, will help to constrain the couplings of the Higgs boson to gauge bosons and 
fermions Q. 

The observation of two forward tagging jets in Higgs production via VBF at the LHC 
is crucial for the suppression of backgrounds ^ 0, ^, |9|, |lO|. In addition to forward jet 
tagging, the veto of any additional jet activity in the central region (central jet veto) leads 
to further suppression of QCD backgrounds such as W~^W~jj , ttjj, and gluon fusion 



Hjj production [g, 11 1. This is due to the fact that the t-channel exchange of quarks or 



gluons tends to radiate harder and more central gluons than in the VBF case. For VBF 
processes, jet activity in the central region is suppressed due to color singlet exchange in the 
t-channel. For the central jet veto (C JV) proposal, events are discarded if any additional jet 
with a transverse momentum above a minimal value, PT,veto, is found between the tagging 
jets [|, H, H, 0, H, 0,0. 

In order to utilize the CJV for the measurement of Higgs couplings, the reduction 
factor, Psurv) caused by the CJV on the observable signal cross section must be precisely 
known. The relevant information is contained in the fraction of VBF Higgs events with at 
least one additional veto jet between the two tagging jets, i.e. we need to know the ratio of 
the 3-jet Higgs cross section, (yjjj, to the inclusive cross section for VBF Higgs production 
with two tagging jets, ajj. The survival probability for the Higgs signal is then given by 
-fsurv = 1 — ^jjj/^jj- Perturbative survival probabilities for the CJV have been calculated 
for the Higgs boson signal and background processes using LO matrix elements j^, 13|. 
The cross section for the VBF process pp — > Hjjj is proportional to Og at LO, which leads 
to substantial theoretical uncertainties (scale variations of 30% or more) . Even though the 
effect on the survival probability is mitigated by the smallness of crjjj/ajj (about 0.1 to 0.2 
for veto thresholds PT,veto ~ 20 GeV), a more reliable prediction requires a calculation of 
the NLO QCD corrections to the Hjjj cross section. We have performed this calculation 
and report on the results in this paper. We do not consider additional reductions of the 
survival probability due to underlying event and pile-up effects. An assessment of these 
effects is best performed after first LHC data have become available. 

A full NLO QCD calculation of the process pp HjjjX involves virtual corrections 
with hexagon diagrams and would be truly challenging. As we explain in Sec. |^, all pen- 
tagon and hexagon contributions are color suppressed by a factor 1/(A^^ — 1) in an SU(N) 
gauge theory, and they are further suppressed by the kinematics of the VBF process. For 
a prediction of the survival probability of the Higgs signal at the few percent level, com- 
mensurate with the knowledge of the VBF cross section for Hjj production at NLO and 
expected experimental accuracies at the LHC, these contributions are completely negligi- 
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ble. We therefore perform the calculation by systematically neglecting the gauge invariant 
subsets of diagrams which involve t-channel gluon exchange and which lead to pentagons 
and hexagons. Similarly, we neglect identical fermion effects for four-quark final states. In 
Sec. ^ we more explicitly specify and justify these approximations and we briefly describe 
the calculation of the LO and NLO matrix elements for Hjjj production. 

Section ^ deals with phenomenological applications of the parton-level Monte Carlo 
program which we have developed. We consider the Hjjj cross section at NLO after typical 
VBF cuts and discuss the reduction of the scale dependence of relevant distributions. We 
show that the scale dependence of Psurv is reduced to about 1% by including the NLO 
QCD corrections to the three jet cross section. Conclusions are given in Sec. Explicit 
formulas for the virtual corrections and for finite collinear terms from initial state radiation 
in gluon initiated processes are collected in two Appendices. 
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2. The NLO Calculation and 
Approximations 



H 



H C H 



The cross section for the leading or- 
der process pp — > Hjjj, via VBF, h 

has been previously calculated as the 
NLO real emission correction to Hjj ^ ' 
production in Refs. ll|,|l^,|o|. The 

relevant Feynman graphs are depicted ° • — 

in Fig. |l|: one needs to consider the 
0{a^as) subprocesses qQ — > qQgH 
and crossed subprocesses with vec- ^ 

tor boson exchange in the t-channel. ^ b 

We explicitly exclude s-channel weak 
boson exchange and thus set aside 

1- , 11 • Trzj- Figure 1: Lowest-order Feynman graphs for _ff 77 7 

higgsstrahlung processes, i.e. V -tlj *' j t, f jjj 

production with subsequent decay 

V — > jj. In the following, higgsstrahlung 

is viewed as a separate process and we also neglect any interference of VBF and hig- 
gsstrahlung (in the case of identical fermion flavors) since these interference effects are 



very small in the phase space region relevant for VBF observation at the LHC [21, 22]. 

In order to clarify our notation and the approximations in our calculation, let us 
start by considering the form of the Born amplitude for the qQ qQgH subprocess. By 
A^3(lg, 2q, 3g, Cg, 6(3) we denote the matrix element for the parton level process 

q{Pa) + QiVb) ^ q{pi) + Q{P2) + 9{P3) + H{P), (2.1) 

shown in Fig. |l]. Two distinct color structures contribute to this Born matrix element with 
three final state colored partons, 

M3{lq,2Q,3g,ag,bQ) = A3(lg,3g,ag;2Q,bQ)5i^i^t'^^-^ 

+ A3i2Q,3g,bQ;lg,ag)6^,iJ^^^^. (2.2) 

Focusing on the gluon emission, each of the amplitudes ^3 can be viewed as a Compton 
scattering amplitude for the process Q{ki) — > Q{k2)g{qi)V{q2), defined by 



MB{k2,qi,q2;eue2) = -egX"^'"^' gsi^{k2) W ^^^^^^ 

I (fc2 + qiY 

+ ^-V'(fci)6i^(ft)e2.(c?2) . (2.3) 

{k2 + qiY J 

Here, —eg"^^^^^ is the left- or righthanded coupling of the quarks to the weak boson, gg 
denotes the strong coupling constant, Pr = |(1 + ^"7^) is the chirality projector, and ei 
and 62 are the polarization vectors of the gluon and of the weak boson, respectively. The 
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role of the polarization vector for the weak boson is taken by a current, h^^ which, for the 
first two diagrams in Fig. |l[ is given by 



h^{Pbn,P2T2) = 5' 



T2TI, V 



e)9HVvg^./'^''Dv[pl,^]Dv[pl^]^P{p2)rPr,i^iPb) 



(2.4) 



with pijf: = Pi — Pj — Pk and pij = Pi — pj, while Dv[q'^] = — My] is the weak boson 
propagator, which, in our calculation, only occurs with space-like momentum. In terms of 



the Compton amplitude of Eq. (2^) the ^3 are then given by 



^3(1,, 3g,ag; 2q, bg) = MB{pi,P3,Pal3; £3, HPbn,P2T2)), 
-43(2q, 3g, bQ] Ig, aq) = MB{p2,P3,Pb23; ^3, HpaTa, pm)) . 



(2.5) 



The gQ — > qqQH subprocess is obtained by crossing the initial state quark q{pa) with 



the final state gluon in Eq. (2.2) and dropping the s-channel graphs which result from 
crossing the diagrams in the second line of Fig. |^. The 3-parton matrix elements A4^ have 
been computed using the helicity amplitude method of Ref . (2^ . 

The real emission corrections to VBF Hjjj production consist of four subprocess 
classes with four final state partons. These classes are (a) qQ qQggH, (b) qQ —f 
qQq'q'H, (c) gQ qqQgH, and (d) gg — > qqQQH. The generalization to the crossed 
processes with q ^ q and/or Q — > Q is straightforward. 



3 3 3 

Figure 2: The dominant virtual QCD corrections. The "blobs" correspond to the sum of all virtual 
corrections to the basic Q QgV Compton amplitude and are given more explicitly in Fig. ^. The 
first diagram and the second pair of diagrams in each line form gauge invariant subsets. 

The above subprocesses lead to soft and collinear singularities when integrated over 
the phase space of the final state partons. We use the Catani-Seymour dipole subtraction 
method to regulate these divergences [ p^ and to cancel them against those originating from 
the virtual corrections. The virtual corrections can be divided into two classes of gauge 
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covariant subsets. The first class (depicted in Fig. ^) are graphs in which the internal gluon 
propagator is attached to a single fermion line and which involve up to box corrections. The 
second class (depicted in Fig. ^) are graphs in which the internal gluon propagates between 
different fermion lines, i.e. they contain a t-channel gluon. These graphs only play a role for 
subprocesses with two initial state quarks or anti-quarks. For gluon initiated processes they 
only contribute to the interference of VBF and higgsstrahlung diagrams, which we neglect. 
The interference of the hexagon and pentagon amplitudes with the Born amplitude is color 
suppressed by a factor do = -/V^ — 1 with respect to the interference of box corrections with 
the Born amplitude. We neglect the contribution of the hexagon and pentagon amplitudes. 
However, in doing so we must also consider the color structure of the real corrections and 
drop contributions which cancel the infrared singularities of the pentagons and hexagons. 



Figure 3: Pentagon and hexagon diagrams for the color structure Si^i^t";^^^. An analogous set 
appears with the external gluon attached to the lower quark line. Note that hexagon graphs with a 
three-gluon- vertex correspond to a color structure which cannot interfere with the Born amplitude. 

As an example for these real emission processes, consider the matrix element for the 
subprocess, 

q{Pa) + Q{Ph) ^ q{pi) + Q{P2) + 9{P3) + 9{P4) + H{P) (2.6) 

depicted in Fig. |^ and denoted by A^4(lq, 2q, 3g, 4^, a^, 6q). M.a has the following color 
decomposition in terms of color subamplitudes, A and 

>f4(lg,2Q,3g,4g,ag,6Q) = (t"3t"^)i,i„(5i2i,^(l,,3g,4g,o,;2Q,6Q) 

Ag, 3g, aq\ 2q, bq) 
+ {t''H''')i,^,5,,i^A{2Q,2.gAg,bQ■^<,^aci) (2.7) 

+ {eH''^)i,^,5^,i^A{2QAg, 3g, hq' 1„ aq) 

+ C1C^(1'?'3.'S;2q,4„&q) 
+ ^nV*S.^(la>43,a,;2Q,33,6Q). 



The A terms correspond to both gluons attached to the same quark line, while B terms 
describe emission of one gluon from each of the two quark lines. Abbreviating these am- 
plitudes by 



Ai = A(lq,2,g,Ag,aq;2Q,bQ) , A2 = ^(Ig, 4g, 3g, a^; 2q, 6q) , 
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■^3 = A{2Q,3g,4:g,bQ;lg,ag) , A4. = A{2q , Ag , 3 g , tQ ] Ig , Ug) , (2.8) 

Bi = B{lg,3g,ag;2Q,ig,bQ) , B2 = B{lq,ig,ag;2Q,3g,bQ), 

the color summed squared matrix element can be written as 

\M^{lg,2Q,3gAg,ag,bQ)\^ = dj,C],{\Aif + |^2p + \Asf + |^4p + \Bi\^ + 

+2x Re [AiA*2 + ^s^I] + 2y Re [(^1 + ^2) • (^3 + ^4)* + BiB^]} (2.9) 

with x = l- Ca/2Cf = -l/(iV^ - 1) and y = l/dc = 1/(A^^ - 1), where the explicit value 
is given for an SU{N) gauge group, with do = -/V^ — 1 and dp = N . The term in Eq. ( p^ ) 
which is proportional to y leads to a soft divergence when integrated over the phase space 
of the soft /collinear parton, which is in fact canceled by the corresponding soft divergent 
hexagon and pentagon graphs. Since we neglect the latter, for consistency, we also need 
to set y = in Eq. ( |2.9| ). The association of the y-terms with the hexagon and pentagon 
diagrams of Fig. ^ is made clear when recognizing that these are all the contributions 
where the same gluon is attached to both an upper and a lower quark line. The y-term and 
the interference of hexagons and pentagons with the Born amplitude are not only color 
suppressed by a factor 1/{N'^ — 1), they are further suppressed because the interfering 
amplitudes are never large simultaneously when typical VBF cuts are applied. Consider, 
for example, the Bi and B2 amplitudes in Fig. ^. Bi is large when and 53 are forward 
(i.e. in the initial qa direction) and (72 and 54 are backward, in the qi, hemisphere. For B2 
to be large, qi and (74 must be forward while q2 and 53 are backwards. These conditions 
cannot be satisfied simultaneously for a large rapidity separation between the highest pT 
jets, which typically will be the two quark jets. The largest interference between Bi and 
B2 and, similarly, between ^1 + A2 and ^3 + Aa is to be expected when both factors in 
the interference terms have similar size, i.e. when both gluons are emitted in the central 
region. For central gluons, however, all contributing amplitudes are suppressed due to the 
gluon radiation pattern of the underlying t-channel weak boson exchange. 

We have estimated the error on the total Hjjj cross section, Acr^^*-^, which we make 
in neglecting the hexagon and pentagon topologies (shown in Fig. ^ and the corresponding 
interference terms (y-terms) in Eq. ( |2.9| ). Consider the dominant phase space region where 
one gluon (say 53) is hard and the second one (54) is soft. The soft emission can be 
factorized as an eikonal factor, while the hard part of Bi and Ai + A2 will be given by 
the upper line of the Born diagram of Fig. |l|, i.e. by ^3(1^, 3^, Og; 2q, 6q). Analogously, 
the hard factor in B2 and ^3 + A4, is ^3(2q, 3g, 6q; 1^, a^), corresponding to hard emission 
from the lower quark line in Fig. ||. Approximately, in the soft region, the y-terms plus 
the corresponding virtual corrections, given by the interference of hexagons and pentagons 
with the Born amplitude, are proportional to the product 
^2 ^2 

£(]#Zl)2Re[^3(lg,33,a,;2Q,6Q)^3(2Q,35,6Q;l„S)*] , (2.10) 

integrated over the 3-parton phase space. In Fig. ^, we compare the absolute value of this 
proxy for the full interference terms (dotted blue curves) with the tree level cross section 
(dashed red curve). Shown is the distribution in rapidity for the veto jet (lowest pT parton) 
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A2. 



B2: 



1 . . '"' , 1 



Figure 4: Feynman graphs for the real emission amplitude ^^4(15, 2q, 3g, 4g, a,, 6q) as described 
in Eq. d^l) 

measured with respect to the center of the rapidity of the two tagging jets. In the right 
panel of Fig. ^, the ratio of the two distributions is shown. As a second estimator for 
the neglected terms, we have calculated the full y-terms for the 4-parton final state and 
soft approximations for the hexagons and pentagons, by keeping the infrared divergent 
C-function terms only, according to the prescription of Ref. ||2^. For both contributions 
the full Catani-Seymour subtraction has been implemented, with dipole terms as listed for 
y = 1 /do in Table The resulting curve for 

da^'-'{HR,fiF)/dyrei 

is shown for ^r = fj.p = AO GeV. The ratio, R, reaches a maximal value of ~ 10~^ 
in the central region between the two tagging jets, in agreement with the result for the 
proxy discussed above. We conclude that the y-terms and the corresponding hexagon and 
pentagon contributions give a relative contribution below one permille everywhere in phase 
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space and can safely be neglected. We note that these interference terms are at the same 
level as the interference between gluon fusion and vector boson fusion for pp — > Hjjj. 
In a complete calculation, not only would the hexagon and pentagon graphs need to be 
calculated, gluon fusion contributions would have to be included as well. 




Yrel Yrel 



Figure 5: The distribution in rapidity of the veto jet measured with respect to the center of rapidity 
of the tagging jets. In the left panel, da/dyrei is shown at LO (dashed red), for the proxy given by 
Eq. ( |2.10| ) (dotted blue), and for the NLO color suppressed contribution in the soft approximation 
(solid green). The right panel depicts the ratio, R, for both the proxy (dotted blue) and the NLO 
color suppressed contribution (solid green). 



In addition to the y-terms discussed above we also neglect any interference terms for 
identical fermions in our simulations. These terms are color suppressed by a factor 1 /N and 
can only contribute when fermion helicities are the same. For charged current contributions 
we have determined the size of these interference terms for 4-quark final states and have 
compared them to the charged current contribution to the LO 3-jet cross section. We 
find a relative contribution of 7.5 • 10~^ within the cuts of Section |^: also these "Pauli 
interference terms" are truly negligible. With these approximations, the f ortran code for 
the real emission matrix element squared was generated with the help of MadGraph |25|. 



The 4-parton phase space integral of the squared real matrix elements suff'ers from 
soft and collinear divergences. The dipole subtraction method of Catani and Seymour 
provides a means to regulate these divergences [^]. In the Catani-Seymour formalism the 
NLO corrections consist of three pieces: (a) the contribution of the dipole subtracted real 
corrections, cr^^^, (b) the contribution of the finite virtual corrections, a^^^ , and (c) a 
piece resulting from the factorization of collinear singularities into the parton distribution 
functions, cf^coY • example, consider the process, qQ qQggH, in the y = case in 

Eq. ( |2.9| ). The subtracted cross section for this process takes the form, 

Cr^^^iqQ ^ qQggH) = [ dXa [ dXbfg/p{Xa, f^F)fQ/p{xb, f^F)-^d<^5{pa,Pb) 

Jo Jo ^'S 

• ||7W4(lg, 2q, 3g,4:g, Ug , ) | ^j^^ (i^l , P2 , P3 , P4 ; Pa +Pb) 
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real subprocess 


dipole factors 


(lg,2Q,3c,,4g,ag,6Q) 
y = U m iijq. (J^.^j 


^14,3) ^13,4, T^34:,l,T^u,T^i3 
T^a 'na4 T^a3 'na4 T^a3 

2^24,3, 2^23,4, 2?34,2, 2^24' ^23 
-T-jft TlW 7^63 'nM -7^63 

■^-^34' -^2 ' ■'^2 1 -^3 ' -^^4 


y = l/dc in Eq. (P) 


pa3,6^ pfe3,a^ pa4,fe^ pb4,a^ 

-^14' ^2^' ^1^' 2^24,1) 21^14,2 
2^23' -^13' -^2 '2^1 1 ■^'23,1, •^'13,2 


Ugi ^Q) -Jg'' % ' ^^9' "Qi 


■t^34,l, ■^'31,2, '^'34, ■'^34 
-nal -n'Jl 71^2 'n62 


(lg)2Q,3q,4g,Og,6g) 


vf,vf,vf,vf 

2^14) 25^4, 1^14,3) 2^34,1 


(lg,2Q,3g,4Q,ag,6g) 





Table 1: Dipole factors for the real emission corrections to Hjjj production. The y = 1/dc line 
gives the additional dipole factors which are needed for the qQ — > qQggH process when the y terms 
in Eq. (|2.9|) are not neglected. 



{pi,..Pij,Pk,--PA]Pa,Pb) (2. 

(3), 



pairs k^i,j 
i,3 



Y [^ij(Pi'P2,P3,Pi,]Pa,Pb)Fy'{pi,..pij, ..,P4]Pa,Pb) + (a ^ 6) 

pairs 

YY [^k'iPl^P2,P3,P3;Pa,Pb)FPipi, ..Pk, ■■■,P4.;Pai,Pb) + (a ^ 6) 



where the 2?ij,fc etc. are the dipole factors as defined in Ref. [26|, is the 5-particle phase 
space measure and s = {pa +Pb)'^ denotes the center-of-mass energy. A complete list of the 
dipole factors in Eq. ( p. 121) is shown in Table ||. Notice, that we do not need to consider 
dipole factors for which there is an initial state singularity with an initial state spectator 
for the case of y = because in this approximation radiative corrections to the upper and 
the lower lines in Fig. |^ effectively decouple. We also show in Table |l| dipole factors for 



quark-gluon and gluon-gluon initiated processes. The functions Fj^^ and Fj*' define the jet 
algorithm for 4-parton and 3-parton final states and must be infrared safe which formally 
^ F^p in any case where the 4-parton and 3-parton configurations are 



means that F 



J 

kinematically degenerate. 

The dipole factors are integrated in d = 4 — 2e space-time dimension over the phase 
space of the soft/coUinear parton. Integrating the dipole factors for the processes, qQ — > 
qQggH and qQ qQq'q'H, lead to the universal singular factor, < 1(e) >. For the parton 
level process 



q{pa) + QiPb) ^ q{pi) + Q{P2) + g{P3) + h{p) 



(2.13) 
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we can split < 1(e) > into two pieces according to, 

< 1(e) >= Cf (Xi(e) +T2(e)) • (2.14) 

Ti(e) is a piece proportional to the Born-level color subamplitude squared, j^3(lq, 3^, o^; 2q, 6q)|^, 
and is 

Ti(e) = |^3(l,,3„a,;2Q,6Q)|2^M)^,^^ (2.15) 



with 



and, 



7, = Icf, 7, = ^Ca - ^TkTV/, (2.16) 



Here Sij = 2pi ■ pj with i = l,2,3,a or b. Tr = 1/2, Ca = N, and Cp = (N^ - l)/(2iV) 
in SU{N) gauge theory. The number of flavors is Nj = 5. 12(e) is obtained from Xi(e) 
by interchanging the quark labels, a b and 1 <-5- 2. The 1/e^ and 1/e divergences cancel 
against the virtual corrections shown in Fig. ||. 

In our approximation there are two distinct color structures that contribute to this 
virtual matrix element, M.™^{lq, 2q, 3^, a^, 6q), with three final state colored partons, 

Mr'{lg,2Q,3g,ag,bQ) = Ar\lg,3g,a,;2Q,bQ)6i,i,t^^^^ (2.18) 

+ Ar\2Q,3g,bQ;l„aq)6,,,j'^X- 

The interference between the virtual and Born three parton amplitudes takes on the fol- 
lowing form upon summing over final state colors and averaging over initial state colors. 



2 Re[Mr'Ml] = Cf (2 Re[Ar\lg, 3„ a,; 2q, bQ)AUU, ^9, 2q, ^q)] (2.19) 



colors 



2 Re[^f *(2q, 3g, 6q; Ig, a,)^^(2Q, 3g, 6q; Ig, a^)]) 



We split the virtual corrections shown in Fig. ^ into two classes: the virtual corrections 
along a quark line with only one weak boson attached and the virtual corrections along a 
quark line with a gluon and a weak boson attached. The former, with only a weak boson 
vertex, are factorizable in terms of the tree- level current defined by Eq. ( |2.4D . For vertex 
corrections to the lower line one has 

K.Mn,P2r.) = h^{p,n,p,r,)CF^^ (^)^ f(i^ (-| - 7 - ^) ' (^-^^^ 
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Figure 6: The one-loop QCD corrections for q qgV. 

and similarly for vertex corrections to the upper line. Here fj-u is the renormalization scale, 
and Sfe2 = 2pfe ■ p2 is the weak boson virtuality for massless quarks. 

The second class of diagrams, corresponding to the "blob" in Fig. |2| and shown explicitly 
in Fig. ^ are the virtual QCD corrections to the Feynman graphs where a gluon g and an 
electroweak boson V ( outgoing momenta qi and ^2) are attached to the same fermion line. 
The kinematics is given by 

Q{ki)^Q{k2)+g{qi) + V{q2), (2.21) 
where kf = k2 = qf = and momentum conservation reads ki = k2 + q2 + Qi- As in 
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Ref. 1 27 1, it is convenient to use the Mandelstam variables for a 2 — > 2 process which is 



taken to be qq gV . The Mandelstam variables are thus defined as 

s = {ki- kif = {qi + g2)^ 

t=iki-qif = ik2 + q2f, (2.22) 
u= {ki -q2f = {k2 + qif. 

The gluon polarization denoted by ei(gi) is transverse, i.e. ei • gi = and this permits 
simplifications in the virtual amplitudes. The electroweak boson V is always virtual in the 
calculation. Its effective polarization vector eg ((72) corresponds to the tree level fermion cur- 
rent /i^. Due to the emission of the Higgs boson off the t-channel vector boson propagator, 
this fermion current is not conserved. Hence, terms with €2-^2 rnust be kept. However, elec- 
troweak gauge invariance of the amplitude is preserved, i.e., A^^^g ~ We have computed 
the virtual amplitude in the conventional dimensional regularization scheme (CDR) and 
used Passarino-Veltman reduction vad = 4 — 2e spacetime dimensions to reduce tensor loop 
integrals into scalar loop integrals |2^. The virtual amplitude My = -Mv(^2, Qi,'l2] ^i, £2) 
for Q{ki) ^ Q{k2)g{qi)V{q2) is 

My = Mb '''[^'^^ A U( + ^^"^^ ' 



Air r(l - e) \ 2 V V y V - W / e 



2Cf VV J \ ~^ J \ ~^ J J ^ ^ 
+ 2 ('^) ^ (-^ - ^) + F[-s, -t, -u) - ^C^ - 8C, \ (2.23) 



+ M 



V , 



where 



Fi-s, -t. -u) = ^ ^1„^ ^-jj + J " 2*^- - '^^^ 

^l(C. - 2Cri in (^) + (hr,N, - \C.) (l„ (^) + 1„ (^)) . (2.24) 

The finite part My = My {k2,qi,q2] (-1,(^2) is given by 

My = '^^{-e)gX^-Q^gs (2.25) 
• I (C^ - \ca){M^^^ ( A:2 , ft , 92; ei , £2) + At^.^) (feg , ft , 92 ; ei , £2) } 
- hjAMf\k2,quq2]ei,e2] 

Results for physical kinematic regions can be obtained through the analytic continuation 
of Eq. ( |2.23| ) by the replacement of the time-like invariant by s — > s -|- iO^ , t — > t -|- iO"*" , 
or n — > n -|- iO"*". The ivr factors which result from the analytic continuation vanish upon 
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interfering the virtual amplitude with the Born amplitude. The analytic continuation for 
any double logarithms is dealt with automatically by the f ortran code for the finite part 
of the virtual amplitude Aiy given by Eq. ( 2.25| ). 

The A4r^ for i = 1, 2, 3 are finite and can be expressed in terms of the finite parts of 
the Passarino-Veltman, Bq, Cq, and Dij functions, which we denote as Bq, Cq, and Dij. 
Analytic expressions for the are given in Appendix 

We can build the virtual color subamplitudes of Eq. ( |2.18| ) out of the two classes of 
virtual corrections discussed above. The virtual color subamplitudes are then 

AY^{lg,3g,ag;2Q,bQ) = Mb{pi,P3, Palsies, Kirt{pbn,P2r2)) 

+ Mv{pi,P3,Pal3;(^3,HPbn,P2T2)) , (2.26) 
^3"'*(2Q,3g,6Q; lg,Og) = MB{p2,P3,Pb23;(^3,Kirt{Para,Pin)) 

+ Mv{p2,P3,Pb23; (^3, h{PaTa,PlTl)). (2.27) 

For the following we adopt the following abbreviations, 

-43,1a = -4.3(lg,3g,ag;2Q,6Q), ^3^26 = -43(2q, 3g, 6q; Ig, a^) (2.28) 
■^3% = ■^T\'^q,^g,aq;2Q,bQ), Al% = -43"*(2q, 3<;, 6q; 1,, a^). 

The color decomposed interference of the Born and virtual subamplitudes is then, 

colors 

1 f f ATTniy ^ f A7rfi%y\ ^ Ca Ig. 



2 VV ^13 / \ Sa3 J J e 



2Cf\\S23J \ Sa3 J \ Sal / / e 

+ 21 [^^X +(^^W 



Sal J \ Sb2 J J ' € 

— Ca - IQCf + -F(Sal, Sa3, S13) 




+ 2 Re[At\lq,3g,aq;2Q,bQ)Al,J , (2.29) 

with Af^{lq,3g,aq;2Q,bQ) = Mv{pi,P3,Pai3;e3,h{PbTb,P2r2))- A similar expression for 
2 Re[^3*2^^3 2b] is obtained by making the replacements, a b and 1 ^ 2, in Eq. ( 2.29 ). 

Summing together the contributions from Eq. (2.15) and Eq. ( 2.29| ) yields the finite 
3-parton NLO cross section 

a^^^iqQ ^ qQgH) = / dxa dxbfq/pixa, fJ-F)fQ/p{xb, fJ-p) 

Jo Jo 

X ^d'^4{Pa,Pb)FP{pi,P2,P3,P;Pa,Pb) (2.30) 



|A^b(1„ 2q, 3g, a„ 6Q)n 1 + ^4z^^born 
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+ |^3(1„ 3^, ag- 2q, bQ)\^^^^^^^F{Sal,Sa3, S13) 



+ \-A3{2Q,3g,bQ; lg,ag)p^^^^^^£^— ^F(sb2,S63,S23) 



27r 



+ Cf (2 Re[Ar\lg,3g,ag;2Q,bQ)A*3{l„3g,ag;2Q,bQ)] 
+ 2 Re[lf *(2q, 3^, bq; 1„ ag)^K2Q, 3^, &q; l^, a?)]) } , 



with 



born 



27r2 
3 



+ y ) Ca - y TfiiV/ + 2Cf (2 - . (2.31) 



The remaining divergent piece of the integral of the dipole factors in Eq. ( |2.12| ) is 
proportional to the P'''^ and P^'^ splitting functions and is factorized into the parton dis- 
tribution functions. The surviving finite collinear terms are given by 

<^3 CO? (IQSH) = I dXa I dXb-^d^4{Pa,Pb)Fj^\pi,P2,P3;Pa,Pb) 

Jo Jo ^■s 



|A^3(lg,2Q,3g;ag,6Q)|2 



pi, a/ 



1 



+ ^CAfq/p{Xa; ^^F){fQ'|p{xb■, f^F, f^n) " f q"/ p{xi,; ^1 F , m)) 

■ l-^3(2Q,3g,6Q; lg,ag)p 

+ \cAifgfpiXa;fJ-F,fJ^R) - f^fpiXa, HF, fJ^R))fq/p{xb, fiF) 

■ \A3{lq,3g,ag;2q,bq)\^}, 



o2,b 



(2.32) 



and similarly for the anti-quark initiated processes. Here the quark functions f^jpix] fJ-F, IJ-r) 
are given by 



fq/pix'^f^F,m) 



2tt 



dz 



+ 



fq/p /^f) - zfg/p{x; hf) Bl'^{z) 



(2.33) 



+ f,/P (f ;/^^) Cii^i^} + ^^fq/p{x;f^F)Dli{x) 



with kernels 



Bi^^iz) = Cf 
Cl^iz) = Cf 
Dl'Kx) = Cf 



Tr[z +{l-z) jln ^ 



l-z 



In 



^l-pZ 

2pjPi{l - z) 7j 1 



+ TR2zil - z) . 



Cil-z 



-{l + z)\n 



2tt 



2pjPi{l - z) 

fl^Z 



l-z 



lnz+ (1 - z) 



Q --^lln(l-x)+ln2(l-x) 



(2.34) 
(2.35) 
(2.36) 
(2.37) 



, 3 2piPj , 2piPj 
+ - In — ^ + 2 ln(l - X) In — ^ 

2 IJ-F fJ-F 
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af^°(Mz) a^^(Mz) Mz 


Mw 


Gp 




0.118 0.130 91.188 GeV 


80.416 GeV 


1.16639 X lO-VGeV^ 


120 GeV 



Table 2: Standard Model input parameters 



where if parton i is a gluon, Cj = Ca and if parton i is a quark or anti-quark, Cj = Cp- 
Likewise, 7i = 7q if parton i is a quark or anti-quark and 7i = if parton i is a gluon. 
The analogous results for gluon initiated processes are given in Appendix]^. 

We have implemented the QCD corrections for pp — > Hjjj into a fully flexible parton- 
level Monte Carlo program. We have checked the dipole subtraction by verifying that 
the dipole subtraction terms and the real emission matrix elements match in the various 
singular regions. The gauge invariance of the virtual matrix elements has been checked 
numerically for random choices of momenta. The finite collinear counter-terms that re- 
main after the factorization of initial-state collinear divergences have been obtained by two 
independent calculations. We have also introduced a cut, a € (0, 1], on the phase space of 
the dipoles as described in Ref. [^9|. We have checked that the integrated cross section is 
independent of this parameter and have used a = 0.3 in our simulations. 

In all subsequent calculations we use the input parameters for defining Standard Model 
(SM) couplings as listed in Table |2[ Other SM couplings are computed using LO elec- 
troweak relations. Cross sections are computed using CTEQ6M parton distributions |30| 
for all NLO results and CTEQ6L1 parton distributions for all leading order cross sections. 
The running of the strong coupling is evaluated at two-loop order, with as{Mz) = 0.118, 
for all NLO results. For LO results, the running of the strong coupling is evaluated at 
one-loop with as{Mz) = 0.130. In order to reconstruct jets from the final-state partons, 
the kp algorithm |3^] as described in Ref. [32| is used, with resolution parameter D = 0.8. 
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3. Predictions for the LHC 



The goal of our calculation is a precise prediction of the LHC cross section for Higgs boson 
production in VBF with three or more jets. The kT algorithm is used to define jets and 
these jets are required to have 

PTj> 20 GeV, |yj|<4.5. (3.1) 

Here yj denotes the rapidity of the (massive) jet momentum which is reconstructed as the 
four- vector sum of massless partons of pseudorapidity |r/| < 5. 

At LO, there are exactly three massless final state partons. At NLO these jets may 
be composed of two partons (recombination effect) or four well-separated partons may be 



encountered, of which at least three satisfy the cuts of Eq. (SJ.) and would give rise to 
either three or four-jet events. As with LHC data, a choice needs to be made for selecting 
the tagging jets in such a multijet situation. Here the "pr-method" is chosen. For a given 
event, the tagging jets are defined as the two jets with the highest transverse momentum 
with 

p^r^f > 30 GeV, \yf^\ < 4.5. (3.2) 

The non-tagging jets by default are jets of lowest transverse momenta. They do not need 
to satisfy the cuts of Eq. (3^) but must satisfy the cuts of Eq. (|3.1|). 



The Higgs boson decay products (generically called "leptons" in the following) are 
required to fall between the two tagging jets in rapidity and they should be well observable. 
While the exact criteria for the Higgs decay products will depend on the channel considered, 
such specific requirements here are substituted by generating isotropic Higgs boson decay 
into two massless "leptons" (which represent t^t~ or 77 final states) and requiring 

Pt£ > 20 GeV , |%| < 2.5 , ARji > 0.6 , (3.3) 

where ARj£ denotes the jet-lepton separation in the rapidity-azimuthal angle plane. In 
addition, the two "leptons" are required to fall between the two tagging jets in rapidity: 

2/?™n + 0.6<%,, <y2„,-0.6. (3.4) 

Note that no reduction due to branching ratios for specific final states has been included 
in the calculation. 

Backgrounds to vector-boson fusion are significantly suppressed by requiring a large 
rapidity separation for the two tagging jets. Tagging jets are required to reside in opposite 
detector hemispheres with 

and to have a large rapidity separation of 

A%, = |yr'-yr'l>4, (3.6) 

sometimes called "rapidity gap cut" . QCD backgrounds for the Higgs signal typically occur 
at small invariant masses, due to a the dominance of gluons at small Feynman x in the 
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incoming protons The QCD backgrounds can be reduced by imposing a lower bound 
on the invariant mass of the tagging jets of 



rujj = \l {p 



tag 1 ^ ^tag 2y ^ 



(3.7) 



The cross section for Higgs production via VBF in association with three jets or more 
(Hjjj), within the cuts of Eqs. (^]|)-(3/7), is shown in Fig. 0. The scale dependence of 
the NLO and LO cross sections is shown for factorization and renormalization scales, fip 
and hr, which are tied to a fixed reference scale /Uq = 40 GeV, 



(3.8) 



The value /xq = 40 GeV was chosen to minimize the scale dependence of the NLO predic- 
tions and at the same time it provides optimal agreement of the LO approximation with 
the NLO result. 
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Figure 7: Scale dependence of the total cross section at LO and NLO within the cuts of Eqs. 



(3.7) for VBF Hjjj production at the LHC. The factorization scale and the renormalization 
scale iJ^n are taken as multiples, i^^o, of the fixed reference scale = 40 GeV. The NLO curves 
are for = fip — (solid red line), fip = and ijlr = ^/xq (dashed green line), and hf = ■Cmo 
and fiR — S^no (dot-dashed blue line ). The dotted black curve shows the scale dependence of the 
LO cross section for /i/? = fiF = ^Mo- 



The LO cross section depends on both the factorization and renormalization scale. 
For IJ,R = fJ^F = CfJ-o with 0.5 < ^ < 2 the scale variation is +2G% to —19% for the 
LO cross section. The large scale variation is primarily due the fact that the LO Hjjj 
production cross section is proportional to a^. This is in contrast to Hjj production in 
VBF, which only depends on the factorization scale at LO. At NLO three choices are 
shown: (a) = = S, (solid red line); (b) = Cf = ^ (dashed green line); (c) = 1, 
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= C (dot-dashed blue line). Allowing for a factor 2 variation in either direction, i.e., 
considering the range 0.5 < ^ < 2, the NLO cross section changes by less than 5% in all 
cases. 




Figure 8: Rapidity separation in Hjjj production within the cuts of Eqs. (3.1)-( |3.5| ) and Eq. (3.7). 
In the left panel, da/dAyjj is shown at LO (dashed green) and NLO (solid red) for — fij^ ~ 
fiQ — 40 GeV. The right-hand panel depicts the K factor (solid green) and scale variations of LO 
(dotted) and NLO (dashed) results for fin = iip = ^/iQ with £^ — 1/2 and 2. 



Our Monte Carlo program allows the analysis of arbitrary infrared and collinear safe 
distributions with NLO QCD accuracy. In order to assess the impact of the NLO corrections 
we compare LO and NLO results by plotting the dynamical K factor 

Kix) = ^"C/^'^ ^^"^ (3-9) 

for our fixed reference scale of //q = 40 GeV. The stability of the results is represented via 
the scale dependence, given by the ratio of cross sections and dubbed "relative change" in 
the following, 

relative change = da^j^n = f^F = ^f^o)/dx 
d(T3{fiR = fiF = fJ-o)/dx 

We plot results for ^ = 1/2 and 2 with fiQ = 40 GeV for NLO and LO distributions. 

The wide separation in rapidity of the tagging jets is a characteristic feature of VBF 
processes. In the left-hand panel of Fig. |8|the distribution da/dAyjj is shown at LO (dashed 
green) and at NLO (solid red) for Hjjj production. Just as in the NLO Hjj case the 
NLO corrections push the peak towards higher values of rapidity separation Ayjj. This 
strengthens the case for the rapidity gap cut of Ayjj > 4. The K factor (solid green) in the 
right-hand-side of Fig. |8| is strongly phase space dependent. The scale variations ^ = 2^^ 
are significantly reduced by the NLO corrections, from ~ 25% at LO to ~ 10% or less at 
NLO in the relevant region 4 < Ayjj < 7. Similar results are found for the transverse 
momentum distribution of the tagging jets as shown in Fig. 
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Figure 9: Transverse momentum distribution for the softer tagging jet in Hjjj production within 
the cuts of Eqs. (|3.l[)-(pj|). The meaning of the curves is the same as in Fig. H. 
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Figure 10: Invariant mass distribution of the two tagging jets in Hjjj production within the cuts 



of Eqs. (3.1)-(|3.5|) and Eq. (|3.6|). The meaning of the curves is the same as in Fig. 



In Fig. 10 the invariant mass distribution of the two tagging jets is shown for a fixed 
value of renormaUzation and factorization scale, /x/; = /Ui? = 40 GeV. The K factor (solid 
green) deviates from unity by 10% or less for this distribution and this scale choice, i.e. the 
LO result provides for an excellent estimate. The ^ = 2^^ scale variations produce changes 
in the LO distribution of about 30%, however (dotted lines). This uncertainty is reduced 
to the 5% level at NLO (dashed curves). 

When contemplating a central jet veto for the VBF signal, the probability for observing 
three (or more) jets in the final state becomes crucial. With the two leading jets defined 
as tagging jets, one would like to know this probability for emitting additional jets as a 
function of tagging jet distributions. It is given by the 3-jet ratio R = (7^/(72, which we 
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Figure 11: The rapidity separation of two tagging jets for nih — 120 GeV within the cuts of 



Eqs. (3.1)- (3. 5) and Eq. (3.7). In the left panel, da/dAyjj is shown at NLO (solid histograms) and 
LO (dashed histograms) for Hjj and Hjjj production with a fixed scale = /i_F = 40 GeV. In 
the right panel, 3-jet ratios, R{Ayjj) are shown at LO (dashed) and at NLO (solid) for pi^j = ^ip = 
20,40, and 80 GeV. 



define for arbitrary distributions as 



For both NLO and LO 3-jet ratios, the distribution for Higgs plus two jet production in 
the denominator is computed to NLO accuracy since this provides the most accurate cross 
section estimate. For these Hjj distributions, the NLO parton-level Monte Carlo program 



described in |19] is used with renormalization scale and factorization scale set to the mass 
of the Higgs boson, rrih. The numerator corresponds to the analogous distribution for Higgs 
plus three jet inclusive events (VBF Hjjj production) for which we explore LO and NLO 
predictions and different scale choices. 

Let us start by considering the scale variations of the 3-jet ratio as a function of the 



rapidity separation of the tagging jets, x = ^yjj, (in Fig. 11) and of the invariant mass, 
X = rrijj, of the two tagging jets (in Fig. |l2|). The left-hand panels show the distributions 
for 2-jet inclusive and 3-jet inclusive events as predicted at LO (dashed histograms) and 
NLO (solid histograms) for a fixed scale ^u/j = fip = 40 GeV. The right-hand panels 
then give the corresponding 3-jet ratios for three choices of scales, = fip = 20,40, and 
80 GeV. 

The 3-jet ratio decreases with increasing rapidity separation of the tagging jets. This is 
largely a kinematic effect: additional radiation in VBF events is mostly emitted outside the 
rapidity range set by the two tagging jets. Thus, the available phase space for additional 
jets diminishes rapidly as Ayjj increases. While typical 3-jet ratios are around 0.2, the 
LO ratio R^'-' {Ayjj) (dashed curves) reaches values up to 0.7 at low values of Ayjj. The 
corresponding NLO ratio is significantly lower, around 0.4. The reason is that at NLO the 
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Figure 12: The dijet invariant mass of two tagging jets for mh — 120 GeV within the cuts of 
Eqs. (3J)-(3^). The meaning of the curves is the same as in Fig. ^ 



separation of the tagging jets increases somewhat. When normahzing to the NLO Hjj 
cross section R^'^ is enhanced in the Ay^j = 3 region, where cross sections are very small 
due to the rrijj > 600 GeV cut. There is no such effect for R^^^. One also notices that 
for higher values of Ayjj the scale dependence decreases, becoming insignificant at NLO 
in the phase space region with typical VBF cuts (Ayjj > 4). 

Similar threshold effects appear in the nijj distributions of Fig. |l^: large scale varia- 
tions at NLO are confined to the low nijj region with negligible cross section due to the 
cuts. The 3-jet ratios decrease somewhat at large values of the dijet invariant mass. How- 
ever, the effect is not as strong as in the Ayjj distribution. Particularly striking is the 
reduction of the scale uncertainty when going from R^^ (~ 30%) to R^^^ (5 to 10%). 

Veto jets are typically defined to be non-tagging jets that reside in the rapidity region 
between the tagging jets. In addition to the cuts of Eqs. ( |3.lD -(^77|), we employ the following 
definition for the veto jets. 



PtT > PT,veto : 



tag 1 tag 2 



(3.12) 



For 4-jet events it is possible to identify two veto jets. In this case, we order the veto jets 
according to their transverse momentum with p^j° ^ > p^j° ^. In the following we take 
PT,veto = 20 GeV unless stated otherwise. 

On the left-hand-side of Fig. 13 the rapidity distribution, da/dt/rei, of the highest pT 



veto jet is shown. Here the rapidity is measured with respect to the average rapidity of the 
tagging jets, 



Vrei = yf ° - {yf' ' + ^ ')/2 . (3.13) 



The two histograms correspond to the LO (dashed green) and NLO (solid red) distributions 
at a scale fiR = fip = GeV. The suppression of jet activity in the center, near Urei = 0, 
is even more pronounced at NLO than at LO, i.e. the higher order corrections strengthen 
the rapidity gap features of VBF events. This is refiected by the K factor (solid line in 
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right-side panel of Fig. 13) which is greater than one for |yj.ei| > 2 and is less than one in the 
central region between the tagging jets. The right-hand-side of the figure also shows the 
scale variations for ^ = 2^^ relative to the ^ = 1 case: the scale dependence is significantly 
reduced at NLO (dashed curves) and remains largest in the regions of small cross section. 
In the vicinity of yrei = the NLO result varies between —20% and +7% down from a LO 
variation of —20% to -|-24%. In the large cross section regions, near y^-d ~ ±2, the scale 
variations at NLO are a few percent only, a drastic improvement from the LO situation. 
This small scale dependence in the large cross section region will be reflected in small QCD 
uncertainties at NLO for jet veto probabilities. 

The effect is clearly visible in Fig. |lj where the transverse momentum distribution for 
the highest px veto jet is shown for fiR = ^ip = 40 GeV at LO (dashed green) and NLO 
(solid red). The scale variations are largest at high but even at a value of 80 GeV 

the NLO results for ^ = 2^^ (dashed curves) deviate from the ^ = 1 case by only —3% 
to +10%. At LO (dotted) these scale variations are —22% to +31%. The K factor (solid 
green) is close to one but decreases monotonically, i.e. at NLO the veto jet becomes slightly 
softer. 
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Figure 13: The distribution in rapidity of the highest pr veto jet with the cuts of Eqs. (3.1)-(3.7) 



and Eq. (3.12), measured with respect to the rapidity average of the tagging jets. In the left panel, 
da/dyrei is shown at LO (dashed green) and NLO (sohd red) for = fiji — 40 GeV. In the 
right-hand panel the K factor (solid green) and scale variations of LO (dotted) and NLO (dashed) 
results are shown for fifj = fip ~ ^/io with ^ ~ 1/2 and 2. 



Fig. shows the effect of the veto cuts defined by Eq. ( p. 12 ) on the tagging jet invariant 
mass distribution. Both LO and NLO 3-jet ratios are reduced compared to Fig. 12 due to 
the restricted rapidity range of Eq. ( p. 12 ) for the veto jets. Fig. |l6| depicts the distribution 
in rapidity separation of the tagging jets with veto cuts. Again, the 3-jet ratios are reduced. 
However, one also finds a significant shape change of the Ay^-j dependence: the fairly steep 
decrease of the 3-jet ratio with increasing Ayjj becomes much less pronounced at NLO. 
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Figure 15: The invariant mass distribution of the two tagging jets for mh = 120 GeV within the 
cuts of Eqs. Eq. and Eq. ( |3.12| ). In the left panel, da/dnijj is shown at NLO (solid) 

and LO (dashed) for Hjj and for Hjjj production at fiR = fip = 40 GeV. In the right panel, 
the corresponding 3-jet ratios, R^'^{mjj) (dashed) and R^^^{mjj) (solid) are shown for the same 
scale choice. 



In Fig. ^ we show the probability for finding a veto jet, 

P.eto = P{pT,.eto) = r dpl^f (3-14) 

as a function of the minimum transverse momentum of the hardest veto jet, pT,veto- The 
scale variations at LO for the absolute veto probability are on the order of up to ±3%. The 
NLO corrections reduce this scale dependence to below the 1% level, i.e. to a negligible 
uncertainty. When imposing a central jet veto, the accepted VBF Higgs production cross 
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Figure 16: Same as Fig. 05^ but for the rapidity separation of the two tagging jets and within the 



cuts of Eqs. Eq. (3.7), and Eq. (3.12) 




PT.veto [GeV] 



Figure 17: Ratio of the 3-jet cross section to the NLO 2-jet cross section for VBF, Pvcto = 
a^ja^^'^ ■ The dashed curves depict LO ratios and solid curves depict NLO ratios for the following 
scale choices: ^ji ~ = 20 GeV (red), = — GeV (green), and ~ = 80 GeV 
(blue). 



section is given by 

<T2(veto) = (1 - Pveto)o-2. (3.15) 

Since, at NLO, Pveto is only about 10% for a veto jet pr threshold of 20 GeV (and lower for 
harder thresholds) the perturbative uncertainty on the SM prediction for the Higgs cross 
section due to a central jet veto is of order 1% only at NLO and hence negligible compared 
to expected statistical errors 
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4. Discussion and Conclusions 



In this paper we have presented the dominant QCD corrections for Higgs production via 
vector-boson fusion in association with three jets. The calculations are implemented in the 
form of a parton-level Monte Carlo program which allows to analyze arbitrary collinear 
and infrared safe distributions with NLO QCD accuracy. 

Our calculation involves several approximations which significantly reduce the com- 
plexity of the virtual corrections. Since we are only interested in phase space regions where 
vector boson fusion processes can be distinguished from QCD backgrounds, we neglect 
contributions which are small once typical VBF cuts, in particular wide tagging jet sepa- 
rations and large tagging jet invariant masses, are imposed. Identical fermion interference 
effects are small after VBF cuts and we have also estimated the contribution from i-channel 
gluon exchange in virtual diagrams (and related real emission diagrams) to be well below 
one permille over the entire phase space relevant for VBF production. Neglecting these 
small contributions, the QCD 1-loop corrections involve only a single quark line and are 
similar in complexity to dijet production in DIS [33|, i.e. they require the calculation of 
box diagrams as the most complex ingredient. 

One reason for the smallness of the t-channel gluon contributions is that they are color 
suppressed, by a factor of 1/{N'^ — 1) in an SU {N) gauge theory, and this feature is generic 
since gluon colors need to be correlated to match the color singlet exchange nature of the 
tree level VBF process. In addition, we find very strong kinematical suppression factors in 
our analysis of the t-channel gluon contributions which can be traced to the characteristic 
gluon radiation pattern in VBF events. It is this kinematical suppression which renders the 
t-channel gluon contributions truly negligible. It would be interesting to find out, whether 
this kinematical suppression persists at higher orders. 

In our phenomenological analysis for the LHC we find that additional jet activity 
between the tagging jets in VBF Higgs production events is even more strongly suppressed 
once NLO QCD corrections are included: K factors go down to 0.7 for jet emission at 
the center between the two tagging jets. This strengthens the case for a central jet veto 
as a background suppression tool. Requiring the absence of any additional jet activity 
of pTj > 20 GeV between the the tagging jets we find veto probabilities for the signal 
of 

-Pveto ~ 10% from this perturbative QCD source. Our NLO QCD predictions for the 
veto probability show small residual scale variations, indicating a relative error on Pveto 
due to higher order effects of 10% or less. This implies that the survival probability 
-fsurv = 1 — -Pveto Can be determined with a perturbative QCD uncertainty of about 1%, 
which is more than sufficient for Higgs coupling determinations at the LHC 0]. 

Beyond the additional jet activity from perturbative QCD radiation, which we have 
analyzed in this paper, additional central jets in VBF events will arise from multiple parton 
scattering (i.e. the underlying event) and from pile-up in high luminosity running. For 
small veto thresholds PT,veto, the contributions from these sources may be as large as 
the perturbative effects which we have considered and need to be estimated independently. 
However, these additional contributions should be independent of the hard scattering event 
and can, hence, be determined from other LHC data, in particular by measuring the jet 
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activity in other VBF processes. What wih be needed on the theoretical side, is a precise 
calculation of the perturbative contribution to the veto probability for these other VBF 
processes, similar to the calculation described in this paper. 
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A. Virtual Corrections 

' — (i) 

In this appendix, we give the expressions for the finite, reduced amplitudes, M.)- that 
appear in Eq. ( 2.23| ) in terms of Bq, Cq, and Dij functions. Here Bq, Cq, and Dij denote 
the finite parts of the Passarino-Veltman Bq, Cq, and Dij functions [28|, and are given 
exphcity below. We write 

M^^Hk2,qi,q2;ei,e2) = V^(A:2){cW(g^i - ^2) + c['Vi + 4'V2 

+ 4'V2(^2+^2Vl}PrV^(fcl), (A.l) 

M^^Hk2, qi,q2; ei, £2) = V^(A;2){cf H^i " ^2) + cf Vi + 4'V2 

+ c^h^Mh+^iy2}PMh), (A.2) 

and, 

M^^Hk2,qi,q2;ei,e2) = V5(A:2){cf Wi - ^2) + cf Vi + 4'V2 

+ cfVVi(^2 +^iV2 + c'f\ty2{h +^2yi}PMh) , (A.3) 

where ei = ei{qi) and €2 = e2(<?2) are the effective polarization vectors for the gluon and 

(7) 

weak boson. The coefficient functions, q , with j = 1,2,3 and i = b,q, 1,2 are given 
below. 



H — 


^^^^(1) 2i?o(t) Tb(gi,t) 
t t 






(A.4) 


— 


Boxf^ +2 €2 • k2T,{qi,t) - 2 €2 ■ 


[Bo{t)- 

'^^ 


Bo{q'2)] 


(A.5) 


^2 — 


Box^^^ + 2 ei • A;iTe(0,i) 






(A.6) 


c(i) = 


BoxW 






(A.7) 


cF = 


^ (2) 2Bo(n) n(g2,n) 
Boxj^ 






(A.8) 


cf) = 


Boxf ^ + 2 £2 • kiT,{qlu) + 2 62 ■ 


M — 




(A.9) 


^2 — 


Boxf'^ +2 ei •A;2Te(0,ii) 






(A.IO) 




Boxf) 






(A.ll) 




= tBoxf ) - 2{tCo{t) + 1) + Bo{t) + ^6(g2^ t) 




(A.12) 


(3)/ \ 


= nBoxf ) - 2(u(7o(u) + 1) + So(n) + n{ql 


u) 


(A.13) 




= Boxf ^ - 2e2 • k2T,{qlt) + 2e2 


mt) - 

■q2 ^_ 


Bom 

ql 


(A.14) 



262 • kiT^{q2,u) - 2e2 • q2 z 2 — 



u-q^ 

= Boxf ^ + -{tCo{t) + l)ei • ki + -(uCo(u) + l)ei • (A.15) 
4=^) = Box^3) (A. 16) 
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The Tj functions are explicitly listed below. 



1 
- ( 

q^B,{q^)}-2q^Co{q\t) (A.17) 



n{q\t) = -—^{2q^[Bo{t) - Boiq^)] + tBo{t) 
t q 



TM\t) = -^{[Bo{t) - 5o(g2)]^i±||! + 2i?o(g') + 1 - 2q^C^{q\t)} (A.18) 



r,(0,t) = i(2Bo(t) + l) (A.19) 

Tfe(0,t) = 5o(t) (A.20) 

Here the coefficients, Box^, for j = 1, 2, 3 and i = b,q, 1,2 are expressed in terms of the 
Passarino-Veltman Dij functions. The Box coefficients with = Dij{qi, ^2, ^2) are listed 
below. 

Boxf ^ = (6^27 + ^Doqi + ^Duqi - Di^qj - 3D23qi + D2iql - 2D25qi 
+ 4526^1 + ^Dot + SDnt - ^Di2t + D2it - D2At + AD25t 

- 4L»26i - ^Dqu - Diiu - -D12U + AD13U + D21U 

- 5D24U + 4:D25u)/2 (A.21) 



Box^^) = -D27ei ■ €2 - 2D3iiei • €2 + 2D3i3ei • €2 - SDriei ■ k2€2 ■ h 

+ SDisei ■ k2e2 ■ ^'2 ~ 4£>22ei • A;2e2 • ^'2 - 8£>24ei • /c2e2 • ^'2 + 12I?26ei • ^262 • ^'2 
~ ~ 3 ~ 3 ~ 

- 4I?36ei • k2e2 ■ /c2 + 4I?38ei • ^262 • k2 + ^^o^^^ ■ 0.2^2 ■ k2 + -D^ei ■ 9262 • ^2 

+ 41)23^1 • 9262 • k2 - 8D25ei • 9262 • k2 + 4D26ei • §261 ■ A;2 - 4Z'3ioei • ^262 • k2 

~ 3 ~ 19 ~ ~ 

+ 4D39ei • ^262 • ^2 - 2-^0^1 ■ ^2^2 " Qx - -^Dx2^\ " A;2e2 ■ q\ + SDiaei • ^262 • q\ 

- 1252461 ■ ^262 • gi + 852561 • ^262 • q\ + 452661 • ^262 ■ ^1 + 453io6i ■ ^262 • 91 

- 453461 • A;262 • gi + 452361 • ^262 • q\ - 452561 • 5262 • q\ - 4535ei • ^262 • gi 

~ 3 ~ 11 ~ ~ 

+ 41)3761 ■ 5262 • gi - 2-^0^1 ■ ^2^2 • g2 - Y-D1261 • A;262 ■ 52 + 4Z)i36i ■ ^262 • 52 

+ 8^)2361 • ^262 • 52 - 41)2461 • ^262 ' 52 - 4Z)266l • ^262 ' 92 " 4D3ioei • A;2e2 • 52 
+ 41)3961 • ^262 • 52 + 41)2361 • 5262 ' 92 " 4^2561 • 5262 • 52 + 4D3361 • 5262 • 52 

- 41)3761 • 5262 • 52 - 2-^12^1 • 6252 + 2^1361 • 625I + --D2361 • 625!) 

- ^52461 • 625^ - 531061 • 625i " -D3361 ■ 625^ + 5376i • 625I 

2 ~ 3 ~ 1 ~ 

+ D3961 • 6252 - -C'061 • 62t - 2-^11^1 ■ + 2^1261 • 62* 

- ^52161 • 62* + ^52461 • 62* + 53io6i ■ 62* - 53561 • 62* 

~ ~ 11 ~ 3 ~ 

+ -D3761 • 62* - -D3961 • 62* - ■ 62lt + -D1261 • 62"^ 
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1 ~ ~ 3 ^ 

- -Diaei ■ e2U + --D2iei • e2U + -D23ei ■ + ^D24€i ■ 

- 252561 • e2U - 52661 ■ e2U - ^aiQei • e2U + 53461 • e2U 

- 53561 ■ 62^ + 53761 • 62^ (A.22) 

Boxf ^ = 2452762 ■ k2 + 20531262 ■ k2 + 2252762 • Ql + 2053ii62 ■ ^1 
+ 12D27e2 ■ 92 + 2OD31362 ■ q2 — 4-D23e2 ■ ^2^2 + 4-026^2 • ^2^1 
+ 253862 • A;2gi - 253962 • k2q2 + 5i262 • qiqj - 5i362 ■ qiqj 

- 352362 • 91 + 52462 • gigi + 252662 ■ qiqj + 253io62 • gigf 

- 253762 • qiq2 - 252362 • 929i + 252662 • 9291 - 253362 • 9292 

~ 3 ~ 3 ~ ~ 

+ 2L»39e2 • 9292 + 2-^0^2 • ^2* + 2-^i2e2 • k2t + 4Di3e2 • k2t 

+ 252562 • k2t + 252662 • ^2* + 2531062 • ^2^ - 253862 • ^2* 

+ 2Doe2 • 91* + 3£'iie2 • qit - -0x262 • qit + 41)1362 • 9it 
+ -02162 • 91* - -D24e2 • 91* + QD25f2 ' 9i* - 2-D26e2 • 91* 

- 2^31062 • 9it + 21)3562 • 91* + 4-Di3e2 • 92* + 4-D23e2 • 92* 

+ 2D2^e2 ■ 92* - 2i:'26e2 • 92* + 21)3762 • 92* - 21)3962 • 92* 
3 ~ 7 ~ ~ ~ 

+ ^Doe2 ■ k2U + -£>1262 • k2U — 2D1362 • k2U + 2D22e2 • k2U 

- 2D2462 • k2U + 2D2562 • k2U — 2D26e2 • k2U + 2D31062 • k2U 

3~ ~ 5~ 

- 20^662 ■ k2U + -Doe2 ■ qiu + -D1162 • 9in + --D1262 • 9i^t 

- 2^)1362 • 9in - 1)2162 • 91 -u + D24e2 • 9i^ - 2D34e2 • 91"" 

~ 3 ~ 3 ~ ~ 

+ 2D3562 • qiu + -£'oe2 • q2U + -Di2e2 • q2U - 2D2362 • 92'" 

+ 2D26e2 • 92"" - 2D31062 • 92'" + 2D3762 • 92"" (A. 23) 



(3) ~ ~ ~ ~ 

Box^ = -I2D2761 • k2 - ^Dsuei ■ k2 - 6D2761 • 92 - 4D31361 • q2 

- SDoei ■ k2ql - 75i26i • ^29! + 25i36i • ^29^ - 25226i • ^29! 

+ 652361 • A;29i - 252461 • A;29| + 45256i • A;29| - 85266i • A;29i 
~ ~ 3 ~ 5 ~ 

- 2D3861 • A;29i + 2-D39ei • ^292 - 2"^°^^ ' ^^^i - ^^u^i ■ 9292 

- 5i36i • 9292 + 352361 • 9292 - -^2461 • 9292 + 252561 • ^292 

^ ^ ^ 3 ,^ 

- 6D2661 • 9292 + 2D33ei • 9292 - 2D3961 • 9292 + ^^^^^ " ^2* 

- 45ii6l • k2t + y5i26i ■ k2t - 252161 • fe* + 252261 • A;2* 



— 6D2561 ■ k2t + 6D2661 • k2t - 2D31061 • k2t + 2D3861 • k2t 

- 2-^oei • 92* - 3Dii6i • 92* + ^D^ei ■ 92* - £>2iei ■ 92* 
+ -D2461 • 92* - 6D2561 • 92* + 6D2661 • 92* - 21)3761 • 92* 
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~ 7 ~ 23 ~ ~ 

+ 2D39ei • q2t + 2-^0^1 ' ^2'" + Y-Di2ei • - GDiaei • A;2U 

— 2D2iei • ^2^ + 21)2261 ■ A;2'u + 10D24ei • — 61)2561 ■ k2U 

— 2D2q€\ ■ k2U — 21)31061 • k2U + 2DsGei ■ k2U + 2Do6i • q2U 

— -D1161 • q2U + 5D1261 • q2U - -D2161 • q2U — 2D2361 • q2U 

+ 3D2461 • q2U + 2D2661 • q2U + 2D31061 • g2'w — 2-D3761 • q2U (A. 24) 

The Box coefficients with D^j = Dij{k2,q2,Qi) are hsted below. 

BoxW = . qi[-8D27 - 8D312 - {Dn - + ^13 - 4^22 + 4524)^1 + Dns - D^s 

+ D13S - AD22S + 4^245 + Dut - Di2t + Di3t - W22t + 4^24*] 

+ 62 • 91 [8^27 + 8^313 + (5ii - 5i2 + 5i3 - 4^22 + ^D2i)ql - Dut - 3Di2t 

+ 3^13* - 4^24* + 4^26*] - 62 • A;2[1653ii - 24^312 - (5ii - 5i2 + ^13 + 4^25 

- 4^26 - 8^310 - 4^32 - 4^34 + 4^35 + 8^36 + 4538)^1 + 5DnS - 5Di2S 

+ Diss + AD21S - AD24S + 4^255 - 4526S - 8D310S + 4535s + 4538S - 4^12* 

+ 45i3i + 4^22^ - 8^24* + 8^25^ - 4^26* - 4^310^ - 4^34* 

+ 4^35* + 4^36*] (A.25) 
Box^^) = ei • A;2[853ii - 2453X3 - (^n + 3^12 - 3^13 - 4^22 + 8^24 - 4^25 + 4^310 

- 41)37 - 4538 + 4539)(j| - 5iiS + 5i2S - 55i3S - 8525.S + 4526S - 4^378 
+ 4539S + 45i2t - 45i3t - 4523* + 4524* - 4525* + 4526* + 453io* - 4537*] 

- 61 • 92 [8527 - 85312 + 245313 + (5ii + 35i2 - 35i3 - 4522 + 8524 - 4525 
+ 45310 - 4537 - 4538 + 4539)9^ + 5ii.S - 5i2S + 55i3S + 8525S - 4526S 
+ 4537s - 4539s - 5ii* - 75i2* + 75i3* + 4523* - 8524* + 4525* - 453io* 

+ 4537*] (A.26) 

Box^-^) = 81)31261 • 62 - 8D31361 • €2 + 8Di2ei • A;262 • ^2 - SDisei • ^262 • ^2 

+ I2D2461 • A;262 • k2 - I2D2561 • ^262 ■ k2 + 41)3461 • ^262 • k2 - 41)3561 • ^262 ■ ^2 

+ D1161 • 9262 • k2 + 7Z)i26i • 9262 • k2 - 7Z)i36i • 9262 • /C2 + 4D226i • 9262 • k2 



+ 8Z)246l 


9262 


/C2 - 41)2561 


• 9262 


• k2 - 8D2661 • 9262 • k2 - 41)31061 • 9262 


■k2 


+ 453661 


9262 


/C2 - D1161 • 


/C262 • 


91 + -C>i26i • /C262 • 91 - -C>l36i • A;262 • 91 




- 452361 


^262 


• 91 - 41)2561 


• A;262 


• 91 + 8D2661 • A;262 • 91 + 4D31061 • A:262 


■ 91 


- 453761 


^262 


91 - 41)2361 


• 9262 


• 91 + 4Z)266i • 9262 • 91 + 4D3861 • 9262 • 


91 


- 453961 


9262 


91 - Z)ii6i • 


A;26i • 


91 + 5D1261 • A:26i • 91 - 5D1361 • A;26i • 9i 


+ 852261 


^261 


91 - 4D2561 


• ^261 


• 91 - 4D2661 • /C261 • 91 - 4D31061 • A;26i 


• 91 


+ 453661 


^261 


91 + 4Z)i26i 


• 9261 


• 91 - 4D1361 • 9261 • 91 + 8D2261 • 9261 • 


91 


- 852661 


9261 


91 + 4D3261 


• 9261 


• 91 - 4Z)386i • 9261 • 91 - 41)31061 • 629! 




- 253261 


629! + 2D3661 • 629! + 2D3761 • 629! + 4-D3861 • 629! 





-32- 



- 2Ds9ei ■ €2ql - 2-^11^1 ' ^2-5 + ^D^ei ■ 628 - --DiaCi • 628 

- 252561 ■ 625 + 252661 ■ 62S + 2^31061 • 62S - 2^3761 • 625 

- 2^3861 ■ 62S + 2^3961 • 62S - 2^2261 • 62* - 2^2361 • 62* 

+ 4^2661 ■ 62* + 4^31061 • 62* - 2^3661 • 62* - 2^3761 • 62* (A.27) 
BOX^^) = -4^27 - 12^312 + 12^313 + 4^31091 + 2^3291 - 2536^1 - 2537^1 - 4538^1 + 2539^1 

- 2Dos - Dns - D12S + D13S + 2525S - 2^265 - 2^3105 + 2^375 + 2538S 

- 2^395 + 2^22* + 2^23* - 4^26* - 4^310* + 2^36* + 2^37* (A.28) 

The Box coefficients witfi Dij = Dij{k2,qi,q2) are listed below. 

Boxf) = -62 • gi[8527 - 8^312 + 24^313 - 4(523 - D26 + 533 " D39)ql + 5iiS - 5i2S 
+ 55i3S + 8525S - 4526S + 4537S - 4539S - 5ii« - 75i2?X + 75i3M + 4523^ 

- 8524W + 4525n - 453ion + 4537U] - 63 • 92 [§527 + 165313 - 4(523 " ^26 + 533 

- 539)^1 + 45i3S + 4523S + 4525S - 4526S + 4537S - 4539S - 5ii7i - 35i2li 

+ 35i3U + 8523M - 4524M - 4526tf - 453iotf + 4537«] - 62 • ^[-85311 + 245313 

- (5ii - 5i2 + 5i3 + 4525 - 4526 + 4533 - 4539)^1 + 5iiS - 5i2S + 55i3S 
+ 8525S - 4526S + 4537s - 4539S - 45i2U + 45i3U + 4523^ - 4524'U + 4525^ 

- 4526n - 453ion + 4537tx] (A.29) 

Boxf = -61 • g2[-8527 - 85313 + (5ii + 35i2 - 35i3 + 4524 " 4526)^] 

- 61 • A;2 [165311 - 245312 + (5ii + 35i2 - 35i3 - 4523 + 4524 + 45310 

- 4537 - 4538 + 4539)^1 + 55iiS - 55i2S + 5i3S + 4521S - 4524S + 4525S 

- 4526S - 85310S + 4535S + 4538S - 45i2U + 45i3« + 4522« - 8524^ + 8525li 

- 4526n - 45310U - 4534U + 4535« + 4536u] (A.30) 

Boxf ) = -4527 - 125312 + 125313 " 2533gi - 2538gi + 4539^1 " 25oS - 5iiS - 5i2S 
+ 5i3S + 2525S - 2526S - 25310S + 2537S + 2538S - 2539S + 2522« + 2D23U 

- Ab2QU - AD^iou + 2D3QU + 2D37U (A.31) 

BOX^^^ = -81)31261 • 62 + 8D31361 • 62 - 81)1261 • A:262 • k2 + 81)1361 • A;262 • k2 

- 12Z)24ei • ^2^2 • k2 + 12D25ei • A;262 • k2 - 41)3461 • ^262 • k2 + 4D3561 • ^262 • k2 

+ -C»ii6i • g2e2 • h - D1261 • 9262 • k2 + D1361 • g2e2 • k2 + 41)2361 • (?262 • k2 

+ 41)2561 • 92^2 ■ k2 - 8D26ei • '3'2e2 • ^2 - 4D31061 • 92^2 • ^2 + 4Z)37ei • '?2e2 • ^2 

- D1161 • /C262 • qi - 7D1261 • A;2e2 • 91 + 7Z)i3ei • k2e2 ■ 9i - 4D2261 • A;262 • 91 

- 8Z)246i • A;2e2 • 9i + 4Z)25ei • k2e2 ■ 9i + 8D26ei • ^262 • 9i + 4Z)3ioei • ^2e2 • 9i 

- 4D3661 • A:2e2 • 91 + 4D23ei • 9262 • 9i - 4D26ei • 92^2 • 9i - 4D38ei • 92^2 • 9i 
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+ ADsgei ■ q2(^2 ■ Qi - ^iiCi • ^262 • 92 - 3-Di2ei • ^262 • q2 + Sl^isei • A:2e2 • q2 
+ 8I?23ei • ^262 • q2 - 4L'24ei • ^2e2 • q2 - 4Z)26ei • k2€2 ■ q2 - 4-D3ioei • ^2e2 • 92 
+ 4L»37ei • /C2e2 • 92 + 4L'23ei • 92e2 • 92 - 4L'26ei • 9262 • 92 + 4D33ei • ^262 • 92 
- 4L>39ei • g2e2 • 92 - 2D33ei • 629! " '^D^^^ei ■ €2qj + 4D39ei • £293 



+ ^Dnei ■ e2S 



^Duei ■ 628 + ^Disei • €28 + 2£'25ei • £2^ 
2D26ei • £25 - 2L>3ioei • €28 + 2L'37ei • £2^ + 2i:)38ei ■ e2S 
2L)39ei • £25 + 2L)22ei • e2U + 2D23ei • 62U - 4L)26ei • £2^ 
4-D3ioei ■ e2U + 2D36ei • e2U + 2L)37ei ■ e2U 



(A.32) 



In the above expressions, the finite Dij functions are obtained by standard Passarino- 



Veltman recursion relations |28| from the finite parts of the basic scalar integrals. For the 
virtual corrections considered, only the one- mass box |34, 35], is needed. Specifically, we 
need the case in which kf = k2 = qf = and g| 7^ 0. Here ki and qi with z = 1, 2 are the 
external four momenta. The one-mass box in the unphysical region, —8 > 0,—t > 0, — > 
is, 



^0(^2,92,91) 



d'^k 



m2 [k^][{k - k^mk -k2- 92)2][(A: - ^2 - 92 - 91)^] 
7r-(/.2)-T(l + 6) 

' t 
2 



2 1 2 1 

8t 8t e 



In 



-9i 



In 



In 



(A.33) 



+ ^o(fc2,92,9i) + 0(e)^ 



where. 



Do{k2,q2,qi) = — 

81 



In 



InM ^ 1 + In' 



+ In 



2 Lio 1 



9i 



(A.34) 



2 Li2 1 



9i 



27r2 
3 



The Mandelstam variables, 8 and t, are defined in Eq. (2.22). 

For the present application, the invariant, is always space-like while the Mandelstam 
invariants, 8 and t, may either be time-like or space-like. Results for physical kinematic 
regions can be obtained by analytic continuation by replacing the time-like invariant by 
t^t + iO+ or 8^8 + iO+. 

In addition, to the one-mass box, we also require expressions for the 3-point and 2-point 
scalar integrals in d = 4 — 2e space-time dimensions. For the 3-point scalar integral. 



CoiPi,P2, (Pi +P2? 



d'^k 



1 



i7r2 [-A:2 - iO+][-{k + pi)2 - iO+] 
1 



(A.35) 



[-{k+pi +P2y-io+y 

two cases are needed. Here pi and p2 represent the external outward flowing four momenta. 
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(a). For the two-mass triangle, either, pf = ov P2 = and p| = (pi + ^2)^ 7^ 0. 



C7o(pf,0,p^) =7r-^(/i^)-T(l + e) 



In 



-pI - Pi 



-pI - iO+ 



In 



Co{plpl) + O{e)] 



Coiplpl) = ^ 



1 1 



In^ 



-pI - iO^ 



In^ 



2-pI-pI \ V 

(b). For the one-mass triangle, p\ = p\ = ^ and p\ = {pi +^2)^ 0- 

1 1 



CoiO,0,pl) =7^~^{^l^r'T{l + e)!^- 



1 , / -p^ - iO+ \ 1 
In ' " ' 



-pI 



- + Co{pl) + O{e)} 



Co{Pi 

The scalar 2-point integral is 
Bo{q^) = 

with 



TT^ 1 11,2 f-pl-iO 



6 -Pi 



+ 



-P3 



2 2 



^2 [_A;2 _ iO+][-(A: - g)2 - iO+] 



vr-(/x2)-T(l + e) 



2 - In- 



^0^ 



B. Cross section formulas 

In this appendix we give cross section formulas for processes of the type, 
giPa) + QiPb) ^ qiPi) + q{P3) + QiP2) + ^^(i^). 



(A.36) 



-pI - «0+ \ \ 1 



2 /^-Ps-^O^ 



(A.37) 



(A.38) 



(A.39) 



(A.40) 



(A.41) 



(B.l) 



Results for the crossed process q{pa) + Q{Pb) q{Pi) + Q{P2) + giPz) + H[P) were al- 
ready given in Section |2|. The finite three parton NLO cross section that results from the 
cancellation of the l/e2 and 1/e poles of the virtual corrections with those of the insertion 
operator, 1(e), is 



O-s^^iaQ ^ qqQH) = / dXa dXbfg/p{Xa,^J-F)fQ/pixb,^J'F) 

Jo Jo 



(B.2) 



1 ro\ 
X ^d^i{pa,Pb)F} '{pi,P2,P3,P;Pa,Pb) 



colors 

+ 2Re[Mr\l„2Q,3g,ag,bQ)M*3{l„2Q,3^,ag,bQ)]] , 
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with 



Mf^lq, 2q, 3q, ttg, bq) = f'^Si^i^MviPl, -Pa,Pal3;ea, h{pbTb,P2T2)) ■ 

The Born level matrix element squared is 

CpN 



\M3{lg,2Q,3q,ag,bQ)f = .^^^^\A3{lq,ag,3g;2Q,bQ)\ 



colors 



with 



Asilq, ag, 3g; 2q, bq) = Mb{pi, 
The finite collinear contribution is 

'^3,coli9Q^QQQH)=J^ dXa dXb^d<^4{Pa,Pb)Flf\pi,P2,P3;Pa,Pb) 
■ {fg/p{Xa;iJ'F)fQ''/piXb;iJ,F,fJ'R) 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



+ 



• j^^-—^\A3{lq,ag,3q;2Q,bQ)y 



with 



In 



+ 



fg/p - zfg/p{xa;nF)\ B'g'^iz) 

Xa)i 



(B.7) 



with kernels, 

Bi'^iz) = Ca 
Cl^^iz) = 2Ca 



l + (l-z)% 2paPi{l-z 



In- 



+ z 



2 j^2p„pi(l-z) 3 1 



1-z 

" 1-z 



- 1 + z(l - z) In 



21- z_ 

2paPi{^ - z) 



1-z 



Inz 



Di'-{x) = 2Ca ln(l - In ^ + 7, In ^ 



+ Ca 



27r^ 50 , 2/-, 
— --+ln^(l-.) 



IJ'F 



(B.8) 
(B.9) 
(B.IO) 

(B.ll) 



+ ^TRNf-h'A-lcAHl-x). 
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